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We investigate the neighbourhood of the chiral phase transition in a lattice Nanibu-Jona-Lasinio model, using 
both Monte Carlo methods and lattice Schwinger-Dyson equations. 



1. INTRODUCTION 

We report the results of an investigation into a 
lattice Nambu-Jona-Lasinio model. For further 
details and references the reader is referred to [1, 

2]. 

The Nambu-Jona-Lasinio model is a purely fer- 
mionic model. The interaction of the particles is 
given by a chirally invariant four fermion interac- 
tion [3] 

- go [ii>(x)^(x)y - {^(x)j,^(x)y]} .(I) 

Our investigation was motivated by the fact that 
four fermion interactions have become popular 
again. They appear for example in technicolour 
models, top-quark condensate models [4] and ex- 
tensions of QED. 

In their original Hartree-Fock calculation [3] 
Nambu and Jona-Lasinio found that for large va- 
lues of the coupling go the fermion mass ij,r is 
generated dynamically and that there exist two 
bound states of a fermion and an antifermion. 
There is scalar state with mass 2ij,r and a mas- 
sless pseudoscalar state which is the Goldstone- 
boson associated with the spontaneous break- 
down of chiral symmetry. 

In order to achieve an understanding of the chi- 
ral phase transition that goes beyond the Hartree- 
Fock approximation we studied the model by the 
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Figure 1. The Schwinger-Dyson equations for the 
fermion propagator. 

Hybrid-Monte Carlo method. Additionally nu- 
merical Schwinger-Dyson techniques were app- 
lied. To enable a direct comparison we studied 
the same lattice version of (1) in both cases. Be- 
cause chiral symmetry is important in the original 
model we used staggered fermions in our lattice 
action, since they allow a continuous chiral sym- 
metry on the lattice. The lattice action reads 
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Figure 2. The Monte Carlo data for the chiral 
condensate and our fit. Symbol shape denotes 
bare mass (diamonds 0.09, triangles 0.04, squa- 
res 0.02 and circles 0.01). 12^* lattices are black 
points, 8^ X 16 white. 

the Monte Carlo treatment the four fermion term 
was rewritten in a quadratic form by introducing 
an auxiliary field 9fj,{x) G [0, 27r): 

- x(^ + //)e-'''''(^)xW] • (3) 

In Fig. 1 we represent graphically the 
Schwinger-Dyson equations for the fermion pro- 
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Figure 3. A comparison between the Monte Carlo 
data for the chiral condensate and the three trun- 
cations of the Schwinger-Dyson equations. The 
lattice size is 8^ x 16, mo = 0.02. 



pagator. The equations as presented there are ex- 
act, but must be truncated to give a numerically 
tractable system of equations. The first approxi- 
mation (known as the gap equation) keeps only 
the term labelled (a). This is the leading term in 
the usual l/N expansion. The next approxima- 
tion is to keep (a) and (b). We refer to this sy- 
stem as the 0(go) equations. Finally we keep all 
three terms but approximate the full four-point 
function by the bare four-point vertex. This ap- 
proximation is referred to as the 0(gl) equations. 

The two methods have different strengths and 
weaknesses. The Monte Carlo yields results wi- 
thout any approximations, but can only be used 
on relatively small lattices. Comparing the re- 
sults of the methods we were able to check that 
the approximations made in the Schwinger-Dyson 
equations are reasonable. Knowing this we could 
confidently use the Schwinger-Dyson technique to 
go to large and even infinite lattices. 

We measured the chiral condensate (xx); 




Figure 4. A plot of the chiral condensate a 
against the renormalised fermion mass. Symbols 
as in Fig. 2, lines as in Fig. 3. 



the renormalised fermion mass ij,r, energies of 
fermion-antifermion composite states and the 
pion decay constant /^ . With these results at 
hand we considered renormalisation group flows. 

2. THE CHIRAL CONDENSATE AND 
CRITICAL COUPLING 

First we looked at the chiral condensate in or- 
der to map out the phase diagram. Using a fit an- 
sate that was based on a modified gap equation we 
determined the critical coupling to be gc Ri 0.278. 
The data and fit are shown in Fig. 2. 

In Fig. 3 we compare Monte Carlo results 
for the chiral condensate with the predictions 
of three different truncations of the Schwinger- 
Dyson equations. The agreement improves as 
more terms are included, and the final curve is 
only slightly shifted from the true results. 

The fermion mass ij,r behaves similarly to the 
chiral condensate. A plot (Fig. 4) of (xx) vs. ij,r 
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Figure 6. The Schwinger-Dyson equations for the 
meson propagator P. 



shows that the data lie on a nearly universal 
curve. The curves are the predictions of the gap 
equation and of the 0(gl) Schwinger-Dyson equa- 
tions. We see that the predictions lie in a very 
narrow band. The differences between the cur- 
ves are caused by finite size effects and a small 
residual dependence on bare mass. 

As can be seen by comparing the white and 
black squares in Fig. 2 finite size effects are rat- 
her severe in the Nambu-Jona-Lasinio model. 
It is therefore fortunate that by analysing the 
Schwinger-Dyson equations we can find a usea- 
ble description of the finite size effects. For large 
lattice size _L, 
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Figure 5. Finite size effects for the chiral con- 
densate on L^ (solid points) and L^ x 2Ls (open 
points) lattices. 
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Figure 7. Energy levels for the tt channel. The 
dotted line is the threshold E = 2//r. 




Figure 8. Energy levels in the tt'/ct channel. Solid 
curves are the tt' levels (the s-wave channel) and 
the dashed curves are the a levels (p-wave). 

(xx)-(xx)(oo) ex {m/Lsf' exp(-//Ri.,) = A .(4) 

This equation can be tested by plotting the chiral 
condensate from differently sized lattices against 
A, as in Fig. 5, (if the equation holds the data 
should lie on straight lines). 



3. SPECTROSCOPY 

In the meson channels we measured propaga- 
tors for the local bilinear operators and for wall 
sources. We compared Monte Carlo results with 
the solution of the Schwinger-Dyson equations 
sketched in Fig. 6. The extra information from 
the wall sources allowed us to find not only the 
ground state energies but also the first excited 
state. In Figs. 7 and 8 we compare our Monte 
Carlo results with the levels predicted (with no 
adjustable parameters) in the Schwinger-Dyson 
approach. 

To understand the meson sector properly we 
need to consider the infinite volume limit. (This 
is particularly important for the levels above the 
threshold E = 2ij,r.) Therefore in Fig. 9 we have 
plotted the tt ground state energies on increasin- 
gly large lattices (8^ x cxd, 12^ x cxd, 20^ x cxd 
and 52^ x cxd). For the largest lattice (thin so- 
lid lines) we have also shown the excited states. 
Also shown are Monte Carlo data from 12'* and 
8^ X 16 lattices. On an infinite lattice the ex- 





Figure 9. Pion energy levels in increasing volu- 
mes. 



Figure 10. The energy levels of the scalar meson. 
(20^ X cxD lattice, mo = 0.02). 




Figure 11. The pion spectral function. In the first 
picture, from the broken phase, the pion is a true 
bound state and the spectral function has a delta 
function at E = m^. In the second diagram the 
pion is a resonance, able to decay into a fermion- 
antifermion pair. 

cited states form a continuum. In the symme- 
tric phase there is a resonance in this continuum, 
which we have plotted as a thick line. On the 
52^ X cxD lattice there is a sequence of avoided 
level crossings These are signatures of the reso- 
nance on a finite lattice. For the pseudoscalar 
meson (tt) we found a light bound state in the 
broken phase. This mass vanishes like ^/mo as 
expected from chiral perturbation theory. In the 
symmetric phase the tt turned out to be a reso- 
nance. The energy levels of the scalar meson (a) 
obtained by the Schwinger-Dyson method are also 
shown in Fig. 10. The cr is a resonance (bold line) 
if go ~ Qc and in the symmetric phase. Deep in 
the broken phase it becomes a very weakly bound 




Figure 12. The spectral function for the scalar 
particle (dashed curve) and the tt' (solid curve). 
The scalar is a narrow resonance near the thres- 
hold. (Note that the two curves are shown at 
different scales.) 

state. In the whole broken phase the a mass is 
about 2iJ,R. 

Using the Schwinger-Dyson equations we can 
also investigate the limit of infinite lattice vo- 
lume, where the closely spaced energy levels lying 
above the threshold go over into a continuum. In 
this case we can define a spectral function p. In 
Figs. 11 and 12 we show spectral functions for the 
pion and scalar channels. We see delta functions 
when a true bound state exists, and resonances 
when a state is heavy enough to decay into two 
fermions. We found the masses and widths of the 
resonances by locating poles on the second Rie- 
mann sheet in complex energy [5] . 

The vector meson (p) also turned out to be 
a resonance. In contrast to the tt and a which 
become massless at the critical point the p does 
not become light. Its mass scales with the cutoff 
when the critical point is approached. 

4. RENORMALISATION GROUP TRA- 
JECTORIES 

Using the results from spectroscopy we ob- 
tained renormalisation group fiows, i.e., contour 
plots of dimensionless ratios of physical quanti- 
ties in the plane of the bare parameters toq and 
go- In particular we looked at the ratios niT^/pR 
and PR/fir, Fig. 13. We found that the lines 
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Figure 13. A comparison between the flows of 
constant m-,; / jjlr (solid lines) and constant jjLrI fir 
(dotted lines) on 8^ x 16 lattices. 
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Figure 14. The m-,; / jjlr flow on an infinite lattice. 
(Dotted line resonance, solid bound state.) 



of constant m-,; / jjlr fiow into the critical point 
while the lines of constant jjLrI fir end on the toq 
axis. The two sets of lines cross everywhere in 
the area around the critical point. Hence there 
cannot be lines of constant physics in this re- 
gion which implies the absence of renormalisa- 
bility. Some caution has to be exercised in in- 
terpreting a fiow diagram measured on a finite 
lattice. As we have seen the ground state energy 
in the pion channel is a true bound state in the 
broken phase but in the symmetric phase it cor- 
responds to an unbound fermion antifermion pair 
with an energy strongly infiuenced by the lattice 
volume. We have therefore repeated the renorma- 
lisation group trajectory on infinite lattices using 
the Schwinger-Dyson equations. The resulting 
fiow diagrams are shown in Figs. 14 and 15. The 
fiow lines cross, so there are no lines of constant 
physics. 
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Figure 15. The jjLrI fir fiow on an infinite lattice. 
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